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The trap environment in which Bose-Einstein condensates 
are generated and/or stored strongly influences the way they 
interact with light. The situation is analogous to cavity QED 
in quantum optics, except that in the present case, one tai- 
lors the matter-wave mode density rather than the density of 
modes of the optical field. Just as in QED, for short times, the 
atoms do not sense the trap and propagate as in free space. 
After times long enough that recoiling atoms can probe the 
trap environment, however, the way condensates and light 
fields are mutually influenced differs significantly from the 
free-space situation. We use as an example the condensate 
collective atomic recoil laser, which is the atomic matter-wave 
analog of the free-electron laser. 



I. INTRODUCTION 

The manipulation of atomic trajectories by optical 
fields forms the basis of atom optics Jj[], a topic of con- 
siderable current interest. Many atom optics experiments 
proceed by simply reversing the roles of light and matter 
from the situation in conventional optics. That is, the 
"optical elements" used to refract, diffract, reflect, focus, 
trap, etc. matter waves are oftentimes — but not al- 
ways — made of light, which transfers momentum to the 
atoms to achieve the desired goal. In many situations, 
it is sufficient to treat the optical fields as constant and 
imposed from the outside. In addition, for low enough 
densities, collisions can be ignored. 

In general, it is necessary to describe both the Maxwell 
field describing light and the Schrodinger field describing 
the atoms as dynamically coupled. Under certain circum- 
stances, one can then formally eliminate the dynamics 
of the atomic field, resulting in effective interactions be- 
tween the light fields. Under a different set of conditions, 
one can alternatively eliminate all or part of the electro- 
magnetic field, resulting in effective atom-atom interac- 
tions (collisions). These are the regimes of nonlinear op- 
tics j| and nonlinear atom optics respectively. But 
they represent limiting cases, where either the atomic or 
the optical field is not dynamically independent, follow- 
ing instead the other field in some adiabatic manner that 
allows for its effective elimination. Outside of these two 
regimes the atomic and optical fields are dynamically in- 
dependent. Neither field is readily eliminated. 

In the past, there have been numerous examples where 
the dynamical coupling between the optical and matter- 



wave fields was important. Much of quantum optics and 
of laser physics deals with such situations. However, in 
the spirit of atom optics, we wish to concentrate only 
on those cases where the center-of-mass motion of the 
material system is central to the problem. An early ex- 
ample of such a situation is the free-electron laser (FEL) 
[^-j?]], where a periodic magnetic field (the wiggler) and 
a running wave optical probe field conspire to spatially 
modulate the density of a relativistic electron beam. This 
density modulation results in an oscillating current which 
amplifies the optical field. As a result, the electron den- 
sity modulation is itself increased, leading to a runaway 
amplification process. A similar mechanism governs the 
collective atomic recoil laser (CARL), || except that in 
that case the relativistic electrons are replaced by atoms, 
and the periodic magnetic field by an optical pump field. 
Here, runaway amplification results from the stimulated 
scattering of the pump field off the density grating im- 
posed on the atomic sample by the combined action of 
the pump and probe fields. Similarly to the free-electron 
laser nomenclature, one would therefore be justified in 
calling the CARL a "free-atom laser." 

Both the FEL and the CARL operate under conditions 
such that the electrons or the atoms can be described as 
classical particles. Indeed, in analogy with the situation 
in optics, one can expect that the wave nature of the 
particles involved is irrelevant if the associated de Broglie 
wavelength is small compared to the characteristic length 
of the "optical elements" involved in their manipulation. 
For the case of the CARL, for example, this implies that 
a classical description of the atomic trajectories is appro- 
priate only for temperatures large compared to the recoil 
temperature. 

The situation becomes in many ways much more in- 
teresting when the wave properties of the atoms become 
essential. Early examples of such situations were dis- 
cussed by Meystre et al, || who studied the diffraction 
of an atom by a quantized field mode, by Haroche et al, 
fllPf who discussed the possible trapping of an atom in 
the vacuum field of a high-Q microwave cavity, and by 
Englert et al, |ll[] who considered the reflection of an ul- 
tracold atom at the entrance of such a cavity. Scully and 
coworkers subsequently extended this work to the the- 
ory of the micromazer, an ultracold-atom version of the 
micromaser. jl2| 

In parallel to these efforts, extensions of the CARL 
theory into the regime of ultracold atoms have also been 



1 



carried out. |13[ They bring to light both the impact 
of atomic diffraction on the threshold behavior of the 
system and the role of the vacuum fluctuations of the 
Schrodinger field on the build-up from noise of the sig- 
nal. Recently, this work was further extended to analyze 
flit the superradiant Rayleigh scattering from a conden- 
sate observed by Ketterle's group. 

These brief remarks touch upon just one of the many 
contributions of Marian Scully yet they clearly illustrate 
his impact on modern quantum optics. It is particularly 
noteworthy that his "ancient" work on the FEL and cav- 
ity QED has now gained a second youth, and is central to 
much of the exciting work taking place in Bose-Einstein 
condensation, atom lasers, and the atom optics of ultra- 
cold atom samples. 

The present paper discusses some new results in 
matter-wave optics that build on the kind of ideas which 
have been central to Marian Scully's research, na mely 
laser physics amplifier theory, and cavity QED [p"7| . 
The new twist here is that instead of tailoring the density 
of modes of the electromagnetic field, we consider a sit- 
uation where it is the matter- wave density of states that 
is manipulated. 

Atomic Bose-Einstein condensates are always formed 
in traps, a feature that has of course important impli- 
cations for many of their properties, such as e.g. the 
dependence of the critical temperature on particle num- 
ber, the quasi-excitation spectrum, etc. In optical exper- 
iments involving the buildup of condensate momentum 
side modes, however, one usually considers a plane-wave 
expansion rather than the actual matter-wave modes in 
the trap. This assumes implicitly that the side modes be- 
have as if in free space. The question, then, is to which 
extent this is appropriate and when a proper description 
of the matter-wave modes of the trap is necessary. 

The problem at hand is similar to the situation in 
cavity QED: for instance, it is known that spontaneous 
emission always first occurs at the free space rate, until 
the wave packet emitted by the atom has a chance to 
probe the cavity boundaries and be reflected back to the 
atom. [p^|Jl9| l This leads to interference effects that can 
result either in enhanced or inhibited spontaneous emis- 
sion, or, for sufficiently high-Q resonators, to reversible 
spontaneous emission. Likewise, it is expected that mat- 
ter waves will behave as in free space for times short 
compared to their propagation time to the trap bound- 
aries. For longer times, however, cavity effects are ex- 
pected to play an important role. In situations where 
matter-wave side modes result from an interaction with 
light, the characteristic atomic velocity is the recoil ve- 
locity v rec — hk/m, which is typically of the order of 
centimeters per second. For trap sizes of 100 microns or 
so, this corresponds to characteristic times of the order 
of tens of milliseconds, after which the atoms will have 
probed their environment and determined that they are 
in a trap. 

This paper illustrates the role of cavity effects in 
matter-wave optics on the example of the ultracold 



CARL. We first briefly present our model, and analyze 
the system's dynamics in the small signal regime in a 
way that exposes the differences between the free-space 
and trap description of the system. We then compare the 
instabilities associated with the two descriptions, show- 
ing the appearance of a new kind of instabilities in the 
"cavity matter-wave optics" analysis. 



II. MODEL 

Our model consists of a Schrodinger field of non- 
interacting bosonic two-level atoms coupled via the 
electric-dipole interaction to two single-mode running 
wave light fields. We consider only the case where the 
detuning between the optical fields and the atomic reso- 
nance is many orders of magnitude larger than the natu- 
ral linewidth of the atomic transition. While single pho- 
ton processes are therefore non-resonant, the atoms may 
still undergo two-photon virtual transitions in which their 
internal state remains unchanged, but due to recoil may 
result in a change in their center-of-mass motion. In the 
far-off resonant regime, the excited state population, and 
therefore spontaneous emission, may be neglected, and 
the ground state atomic field then evolves coherently un- 
der the effective Hamiltonian 
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where m is the atomic mass, V(r) is the trap potential, 
gi and g 2 are the probe and pump coupling coefficients, 
respectively, and K = ki — h 2 is the difference between 
the probe and pump wavevectors. The operator a\ is the 
photon annihilation operator for the probe mode, taken 
in the frame rotating at the pump frequency lo 2i hence 
the energy of the probe mode is the difference between 
the probe frequency u>i and that of the pump. The pump 
is treated classically, and assumed to remain undepleted. 
Thus a 2 is simply a constant, related to the pump inten- 
sity I 2 by | <?2 1 2 |o-2 1 2 = d 2 I 2 /2h 2 eoc, where d is the mag- 
nitude of the atomic dipole moment. We remark that 
we have neglected terms corresponding to the spatially 
independent light shift potential. These terms do not in- 
fluence the dynamics, as all atoms inside the light field 
undergo the same phase shift. Note that this shift mod- 
ifies the dispersion relation for the light fields, hence by 
ignoring it we are taking the index of refraction of the 
atomic sample to be the same as that of the vacuum. 

We assume that the atomic field is initially a Bose- 
Einstein condensate with mean number of condensed 
atoms N, and that this condensate is well described by a 
number state so that the initial state of the atomic field 
may be taken as 
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where |0) is the vacuum state, and 
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is the creation operator for atoms in the condensate state 
ipo (r) , which in the absence of atom-atom interactions is 
the ground state of the potential V(r). 

The Heisenberg equation of motion for the condensate 
field operator is readily derived from (|l|), giving 

^-c = —iuuqCq — i-j= [ d 3 r<y9o(r) [a^e~ ,K ' r 
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where Hloq is the ground state energy, \ — 
|ffi||52||a2|V^V/|A| is the effective coupling constant be- 
tween the condensate and the probe field, and a — 
(5i92 a 2^/l5iHff2||a2||A|)ai is simply the probe annihi- 
lation operator times a phase factor related to the phase 
of the pump laser and the sign of the detuning. 



A. Free-propagation regime 



B. Cavity regime 

The free propagation model is valid for times short 
enough that atoms at the recoil velocity propagate only 
over distances which are short compared to the dimen- 
sions of the initial condensate. As we have shown, in this 
regime the atomic field can be expanded onto momentum 
side modes, which are simply plane waves with a slowly 
varying spatial envelope. A second well-defined regime 
occurs at much longer time scales, when the recoil veloc- 
ity atoms propagate over distances which are larger than 
the trap dimensions. In this 'cavity atom optics' regime, 
the atomic field operator is best expressed in terms of the 
trap eigenmodes {ip n (r)} according to 
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where c„ is the annihilation operator for atoms in mode 
n. With this expansion, the equation of motion for the 
condensate field operator becomes 



From Eq. (|4j) we see that the condensate mode is opti- 
cally coupled to two new states whose field operators are 
given by 
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These states retain the spatial probability distribution of 
the condensate state, but propagate at the recoil velocity 
v r = h[K.\/m, hence we refer to them as 'momentum side 
modes' of the original condensate. With these definitions, 
Eq. (0) becomes 
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Due to the fact that the recoil velocity is extremely slow 
(~ 1 cm/s), the propagation of the side mode wavepack- 
ets can be neglected for reasonably long times, in which 
case the side mode field operators obey 
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where io r — Ti\K\ 2 /2m is the recoil frequency. Wc 
note that we have neglected coupling to higher order 
side modes, an approximation valid in the small-signal 
regime. Lastly, we note that in this free propagation 
model the probe field operator satisfies the equation of 
motion 
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where 5 = lj\ — u>2 is the pump-probe detuning, and we 
again neglect higher-order momentum side modes. 
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where 
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is the matrix element for the optical transition, and we 
have assumed without loss of generality that A nn > is a 
real number. 

We assume for simplicity that the matrix elements 
A nn i is sharply peaked, so that |n — n'\ — m where 
m is some number. This could be the case, e.g. in a 
Fabry-Perot-type matter-wave resonator, where the ab- 
solute value of the momentum is relatively well defined 
for each trap energy level. This allows us to neglect all 
excited states except n = m, which gives 
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In this simplified model, the condensate mode is coupled 
to a single excited trap level, whose equation of motion 
is given by 



dt 
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where fiuj m is the energy of the the mth trap level, and 
we again neglect higher-order couplings. To complete 
this model, we also require the equation for the probe 
field operator, 
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Equations jT^-p^|) differ from ([||9|) in that the two run- 
ning wavepacket momentum side modes have now been 
combined into a single standing wave atomic field mode. 
In the free-propagation regime, which state the atom 
jumps into depends on whether a probe photon was ab- 
sorbed or emitted, a result of momentum conservation. 
In contrast, in the cavity regime both processes transfer 
atoms to the same quantum state, due to the uncertainty 
in the momentum direction of the trap levels. We will 
see shortly that this distinction can lead to remarkably 
different physical effects. A similar situation is famil- 
iar in conventional atom optics, where the diffraction of 
an atom by a standing wave is in general different from 
the diffraction by two countcrpropagating running waves. 



C. Linear response 

The condensate mode is assumed to be initially occu- 
pied by a very large number N of atoms, in which case 
it can be treated classically. Introducing a spontaneous 
symmetry breaking ansatz and replacing the operator fi 
with the c-number order parameter c(t), and considering 
furthermore times short enough that the side mode pop- 
ulations are small compared to N, so that we can neglect 
condensate depletion, we have 
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It is natural to choose the zero of energy to coincide with 
the ground state energy, so that luq — without loss of 
generality. Substituting expression ( |l6| ) into Eqs. (^]|) 
and linearizing in the side mode and probe field operators 
yields a closed 3x3 system of equations, 
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which describes the linear response of the system in the 
free propagation regime. We note that while there are 
six independent operators cL, fi+, d, and their Hermitian 
conjugates, the linearized equations have decoupled into 
two sets of three coupled equations due to the momentum 
selection rules. 

In the cavity regime, there are only four coupled opera- 
tors, c m , c m , a and a}. The equations of motion, however, 
do not decouple and the dynamics is determined by the 
4x4 system of equations 





Cm 




( 









-Xm 




fit 

'-m 


= i 









Xm 


Xm 




a 




\ 


Xm 


~Xm 


-<5 





V 


at J 




Xm 


Xm 





6 



d_ 

dt 



where Xm = xMr, 



( Cm \ 



a 

at 



(18) 



Equation ((F^) for the free propagation regime, and 
Eq. ( f!8|) for the cavity regime can be derived from the 
effective Hamiltonians 
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respectively. Terms like a'& correspond to the genera- 
tion of correlated atom-photon pairs. This is analogous 
to the optical parametric amplifier (OPA), except that 
in that case it is correlated photon pairs that are gener- 
ated. The OPA is currently the primary device for the 
generation of entangled and/or nonclassical states, and 
plays a role in many fundamental experiments in quan- 
tum physics, such as tests of Bell's inequality, quantum 
cryptography, and quantum teleportation. This analogy, 
and the possibility of generating entanglement between 
atomic and optical fields is discussed in detail in pl| , pT| . 

Both linear systems can be solved analytically, all that 
is required are the eigenvalues and eigenvectors of the 
3x3 and 4x4 matrices which appear on the right-hand- 
side of equations dl7| ) and (|l^) respectively. The char- 
acteristic frequencies in the free space regime satisfy the 
cubic equation 
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This equation either has three real solutions, in which 
case the system is stable and experiences only small os- 
cillations about the initial state, or it has one real and 
a pair of complex conjugate solutions. In the latter case 
the eigenvalue with the negative imaginary part corre- 
sponds to an exponentially growing solution, and hence 
the system is unstable. In this case any small input, 
even quantum fluctuations, results in a large output, i.e. 
the system acts as an amplifier for the optical probe and 
atomic side mode fields. 

In the cavity regime, the characteristic equation is a 
quartic equation, given by 
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which is simply a quadratic equation in oj 2 . The solutions 
to Eq. ( p2[ ) fall into three categories. The first category is 
of course when all solutions are purely real, in which case 
the system is stable. A second possibility is that there 
are two purely real and two purely imaginary solutions of 
the form {uj\ =VL,u)2 — —17,^3 = iT.uj^ — — iT}, where 
f2 and T are both real quantities. In this case there is 
only one exponentially growing solution, at the imaginary 
frequency LU4. A final possibility is that the solutions 
are complex numbers of the form {oj\ = Q + iT,uj2 = 
— fi + iT, 0J3 = il — iT, 0J4 = —CI — iT}. This last case is 
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interesting in that there are two exponentially growing 
solutions, U3 and U4, which grow at the same rate T, but 
rotate at equal and opposite frequencies ±0. This leads 
to the temporal modulation of the exponential growth 
rate. 

In Fig. 1 we plot the domains of parameter space cor- 
responding to the various types of solutions. The shaded 
areas correspond to the unstable regions where exponen- 
tial growth occurs. Figure la shows the free propaga- 
tion model, where the stable and unstable regions are 
depicted in the S-\ 2 plane. Both S and x are taken in 
units of uj r . In Fig lb, we plot the corresponding results 
for the cavity model in the 5-\m plane. Here S and Xm 
are taken in units of uj m . In the figures, roman numeral 
I corresponds to a single exponentially growing solution 
and roman numeral II corresponds to two exponentially 
growing solutions. 



From expression (p3|), together with the specified ini- 
tial condition, we find that the mean intensities I\ = 
(cLc_), I2 = (c\c+), and I3 = (a^a) take the form 
Ij(t) = Ij exp(2rt), where 

^■HM' + M'M 2 . (26) 

This expression contains a spontaneous contribution, 
present even in the case a = 0, as well as a stimulated 
term proportional to the injected signal strength |cc| 2 . 
The former is due to the amplification of vacuum fluctu- 
ations in the atomic density and the probe electromag- 
netic field, whereas the latter is due to amplification of 
the injected probe field. 

B. Cavity atom optics regime 



III. MEAN INTENSITIES 

In this section we use the linearized description of sec- 
tion II to derive expressions for the mean intensities of 
the probe and side mode fields, and compare the results 
for the two models. We focus in particular on the 'ex- 
ponential growth regime', which occurs for times long 
enough that all but the exponentially growing solutions 
of section II can be safely neglected, yet short enough 
that condensate depletion can be ignored. We assume 
for concreteness that the atomic side modes begin in the 
vacuum state, while the probe field is initially in a co- 
herent state a, corresponding to the injection of a weak 
laser field into the ring cavity. 



A. Free propagation regime 

In the free propagation regime, the solutions to 
Eq. ( |l7| ) is given by 





= A(t) 



The 3x3 matrix A(t) is given by 

A(t) = Uexp(£W£)U-\ 
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where Uit is the «th component of the fcth eigenvector 
of the matrix on the right-hand-side of Eq. (|l7|), and W 
is the corresponding diagonal matrix of eigenvalues. For 
times long compared to the exponential doubling rate, 
we can neglect all but the exponentially growing terms, 
in which case we have 
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where a%j is a constant which depend only on the control 
parameters 5 and \- 



by 



In the cavity regime, the solution to Eq. ( |l8| ) is given 
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where the 4x4 matrix B is defined in the same manner 
as A, but the eigenvalues and eigenvectors are those of 
the 4x4 matrix on the right-hand-side of Eq. (jlj) . 

Again for moderately long times we keep only exponen- 
tially growing terms, which in region I of Fig. lb allows 
us to take 
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where fry is a constant which depends only on the pa- 
rameters (5 and Xm- With this approximation, we find 



that the mean intensities, 
are all of the form Ij (t) = 
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Ij exp(2rt), where 
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As in the free propagation model, there is a spontaneous 
component, present even for a — 0, and a stimulated 
component, which is present only for a^O. One impor- 
tant difference, however, is that in contrast to the free 
propagation regime of Eq. (p6|), where the intensities de- 
pend only on \a\, in the cavity situation the intensities 
are also sensitive to the phase of the initial probe field. 
We remark that the intensities in region I of the cavity 
model and those in the free propagation model are very 
similar in that the time-dependence is a pure exponen- 
tial. This is no longer the case in unstable region II of 
the cavity model. 

In region II of Fig. lb, there are two independent ex- 
ponentially growing solutions. Keeping only these terms 
gives 
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where Cjj and d%j are again constants which depend only 
on the parameters S and Xm- Due to the counterrotat- 
ing terms, one immediately expects to see oscillations, or 
beating, superimposed onto the exponential growth. 

The mean intensities in region II take the form Ij (t) = 
Ij(t)exp(2Tt), where 

Ij(t) = Xj + 2Re [x 3 e l2nt ] + \a\ 2 (Yj + 2Re [%e i2n *]) 
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Zj = <1, ■■ 

From Eq. (|3l"| ) we see that in addition to purely expo- 
nential terms, there are now contributions with superim- 
posed oscillations. These beats occur even when quan- 
tum noise alone is amplified, as can be seen in the first 
and second terms in Eq. (31). These terms are sponta- 
neous contributions, present even in the absence of an 
injected probe field (a = 0). In contrast, the third term 
is a stimulated contribution which is independent of the 
phase of the injected field and the fourth term, which 
gives additional oscillatory and non-oscillatory terms, is 
sensitive to the phase of the injected field a. 

We remark that when the transcendental equation 
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appearing in that last contribution is satisfied, say at 
a time t = t$, then the phase sensitivity of the so- 
lution disappears. This property recurs for the times 
t = t + mr/fl, where n is any integer. At these times, 
the intensities corresponding to initial states having the 
same \a\ but different phases intersect. 

The growth of the matter-wave intensity is illustrated 
in Fig. 2 for two examples. Figure 2a shows In I/3I versus 
t for the case a = 1 for Xm = w„ and S = — w~ , a 
set of system parameters in region I. Here, there is only 
one exponentially growing solution and no oscillations. 



Figure 2b shows the result for Xr, 



and 5 



This example demonstrates the beating which occurs in 
region II as a result of the interplay between the two 
counterrotating exponential terms. In both figures t is 
taken to be in units of u> m . 



IV. SUMMARY AND CONCLUSION 

The fact that Bose-Einstein condensates of low den- 
sity atomic vapors are generated in traps has impor- 
tant consequences on their intrinsic properties. For in- 
stance, in a harmonic trap the dependence of the num- 
ber of condensate atoms on the critical temperature 
scales as 1 — (T/T c ) 3 , rather than the free-space result 
1 — (T/T c ) 3 / 2 . The condensate ground state and the spec- 
trum of quasi-excitation also depends drastically on the 
presence and form of a trap potential, so much so that 
it is even possible to create stable small trapped conden- 
sates of atoms subject to an attractive interaction. 

In this paper, we have discussed how the trap environ- 
ment also strongly influences the way it interacts with 
light. We have shown that after times long enough that 
the condensate atoms can probe their environment, the 
way condensate side modes and light fields are amplified 
differs significantly from the free-space situation. The 
situation is analogous to cavity QED in quantum optics, 
except that in the present case, one tailors the matter- 
wave mode density rather than the density of modes of 
the optical field. 

We have considered the example of the ultracold collec- 
tive atom laser, or "free-atom laser," which is an atomic 
matter-wave analog of the free-electron laser. We have 
demonstrated how the stimulated scattering of light by 
matter-wave fields leading to their mutual amplification 
is drastically changed in matter- wave resonators, as com- 
pared to the free-space situation. This paper has con- 
centrated on the growth rate of the intensity of the side 
modes of both the condensate and the optical field. Fur- 
ther work will investigate in which way such trap environ- 
ments also modify the quantum statistics of these fields 
and their entanglement. 
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FIG. 1. Exponential instability regions for the free propa- 
gation and cavity regimes are shown in Figs, la and lb re- 
spectively. Region I corresponds to one exponentially growing 
solution, whereas in region II there are two counterrotating 
exponential solutions. In Fig. la, S and \ are taken in units 
of u) r , while in Fig. lb <5 and x™ are m units of uj m . 



FIG. 2. Figure 2a shows In | ^3 1 versus t for the case a — 1 
for Xm = ^>m and S — —uj^f, a set of system parameters 
in region I. Figure 2b shows the result for Xm = '-oTa an d 
S = w„ . This example demonstrates the beating which oc- 
curs in region II as a result of the interplay between the two 
counterrotating exponential terms. In both figures t is taken 
to be in units of ijj m . 
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